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ABSTRACT 



Aims. We explore the sensitivity and completeness of long baseline interferometric observations for detecting unknown, faint com- 
panions around bright unresolved stars. 

Methods. We derive a linear expression for the closure phase signature of a faint companion in the high contrast regime (< 0.1), and 
provide a quantitative estimation of the detection efficiency for the currently offered four-telescope configurations at the Very Large 
Telescope Interferometer. The results are compared to the performances provided by linear and Y-shaped interferometric configura- 
tions in order to identify the ideal array. 

Results. We find that all configurations have a similar efficiency in discovering companions wider than 10 mas. Assuming a closure 
phase accuracy of 0.25 deg, that is typical of state-of-the-art instruments, we predict a median dynamic range of up to six magnitudes 
when stacking observations obtained at five different hour angles. 

Conclusions. Surveying bright stars to search for faint companions can be considered as an ideal filler programme for modern inter- 
ferometric facilities because that places few constraints on the choice of the interferometric configuration. 

Key words. Techniques: interferometric - Binaries: close - Stars: low mass, brown dwarfs - Planetary systems 
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1. Introduction 

The most successful method for detecting faint companions 
around nearby stars is u ndoubtedly the radial velocity (RV) tech- 
nique dPerrvmaril2000l with more than 500 extrasolar planets 
detected up to date. Nevertheless, RV has three major draw- 
backs. First, it cannot be applied to all kind of stars: pulsating 
stars, as well as young, rapidly rotating and/or active stars have 
an intrinsic radial ve locity jitter that gener ally precludes planet- 
search programmes (lUdrv & Santosll2007l) . The discovery of ex- 
trasolar planets orbiting high-mass main-sequence stars (O, B, 
and A spectral types) and young solar-type stars (T Tauri stars) 
therefore remains rare in the literature. Second, RV can only ac- 
cess companions with periods shorter than a few tens of years, 
owing to the limited time span of the observations obtained with 
high-precision spectrographs. Third, RV is inherently insensi- 
tive to the orbital inclination i. This implies that additional astro- 
metric observations are mandatory to unveil the real companion 
mass from the measured quantity M p sin i. For these reasons, the 
development of direct imaging techniques for low-mass com- 
panions is curr ently one of the highest priorities in instrumen- 
tal as tronomy ( Qpp enheimer & Hinkle\i l2009: Absi l & Mawetl 
120101) . 

Direct imaging techniques should be tailored to the typi- 
cal angular separations for which low-mass companions are to 
be found. For instance, the search for young planets in nearby 
star-forming regions (> lOOpc) should be optimised for typ- 
ical linear separations from 0.1 to a few 10 AU for Sun-like 
stars, where planets are supposed to be forming and migrating 
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in their early history (iBodenheimer & L in 2002). This translates 
into typical angular separations of 1-100 milli-arcsecs (mas). 
In a more general context, more than 90% of the sub-stellar 
companions detected to date by all techniques have expected 
appar ent separations smaller than 100 mas (ISchneider et al.l 
2011). This domain is hardly accessible to single-pupil imag- 
ing techniques, even with adaptive-optics (AO) assisted facilities 
dBurrowsl20()5). Nowadays, with 10-m class telescopes, only in- 
terferometric aperture-masking observations routinely break the 
100 mas limit. However, their search space is still restricted to 
angular separations la rger than about 40 mas for dynam ic ranges 
of the order of 500: 1 dKraus et al.ll2"008t lLacour et al.ll201 ll) . 

Thanks to its higher angular resolution, optical long baseline 
interferometry is the ideal tool for exploring separations in the 
range 1-50 mas. The simplest strategy for faint companion de- 
tection with optical/infrared interferometry is to obtain closure 
phase measurements that are directly sensitive to the asymme- 
try in the brightness distribution of the target source and hence 
to possible off-axis companions. However, owing to the sparse 
structure of the point spread function associated with the diluted 
aperture of an interferometer, the depth to which a companion 
can be detected strongly depends on the relative orientation of 
the companion and the interferometric baselines. Consequently, 
the sensitivity limit should be defined as a two-dimensional map 
and/or for various completeness levels. 

The goal of this paper is to provide quantitative estimations 
of the detection efficiency versus the companion contrast and 
separation considering realistic observations. In Sect. |2] we de- 
rive a linear expression for the closure phase signature of a faint 
companion. We also define a simple, general expression for the 
sensitivity limit in terms of companion contrast, which we use in 
Sect.[3]to compute the efficiency of various interferometric con- 
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figurations. We compare the efficiency of the four configurations 
currently offered at the Very Large Telescope Interferometer. To 
add generality, we also discuss the efficiency of three "standard" 
interferometric arrays (two linear and one Y-shaped). In Sect.|H 
we discuss the typical achievable dynamic range for the closure 
phase accuracy provided by existing instruments. 



2. Detecting faint companions with closure phase 

Optical long baseline interferometers provide several observ- 
ables, which are all sensitive to the presence of a faint compan- 
ion: visibilities, differential phases, and closure phases. The lat- 
ter has the main advantage of being, to first order, uncorruptedby 
telescope-specific phase errors, including pointing errors, atmo- 
spheric piston, and longitudinal dispersion due to air and water 
vapour. The noise floor is only limited by the accuracy to which 
the instrumental systematics can be mastered. As a consequence, 
the high-precision closure phase is the favoured observational 
strategy to directly detect faint companions with modern inter- 
ferom eters dVannieretal.1 120061 IZhao et al.ll200l iRenard et al] 
2008). 



2. 1 . Closure-phase signal of a faint companion 

Measuring a closure phase requires the use of an interferometric 
array composed of (at least) three telescopes. Each pair of tele- 
scopes defines a geometrical vector referred to as the interfer- 
ometric baseline. The closure phase is then defined as the sum 
of the phases measured on the three baselines, or equivalently, 
as the argument of the bispectrum, formed through the triple 
prod uct of the mea sured complex visibilities around the trian- 
gle (iMonn ier 2003). Assuming that the stellar diameters are not 
resolved at the considered interferometric baseline lengths, the 
closure phase signature of a binary object can be expressed as 
a function of the three baseline vectors projected onto the sky 
{B\2,B23,Bi\}, the wavelength of observation A, the binary flux 
ratio p (referred to as the contrast in the following), and the ap- 
parent binary separation vector A 



= arg 



where 



<(l + pe' a ^-)(l +pe"**)(l + pe ,a ^) 
\ (1+P) 3 



(1) 



Bi2 A „ B 2 3 A „ B M A 
a n = 2n — ; a 23 = 2n — ; an = 2n — . (2) 

AAA 

The geometrical terms an, «23, and an describe the relative 
orientation of the companion compared to the spatial frequencies 
explored by the interferometric array. We note that an - -(an+ 
023) since the vectorial sum of the three interferometric baselines 
is zero by definition. 

To model the signal of a high contrast binary, we can ei- 
ther use Eq.Q~]with small (<k1) or large (»1) values of p. In the 
first case, the interferometer points toward the primary, while 
in the second case it points toward the companion. Since the 
closure-phase is independent of the pointing position, these two 
approaches are formally identical. Assuming p <k 1, we can de- 
velop and simplify Eq.[T]to be 



; p (sin an + sin 0-23 - 801(0-12 + #23)) ■ 



(3) 



The quantity m = sinai2 + sino'23-sin(a'i2 + Q'23) is referred to as 
the magnification factor because a high value of m corresponds 
to a stronger companion signature in the closure phase signal. 



The validity range of our approximate formula can be ex- 
plored numerically by computing the difference between Eq. Q] 
and[3]for various values of {p,an,a23}- We found that Eq.[3]stays 
within a factor < 1 .5 of the value given by Eq. Q] as long as 
p < 10 , which we consider to be the validity range of our 
work. Therefore, the following results apply only to relatively 
faint companions, and cannot be immediately transposed to bi- 
naries with flux ratios close to unity. 

Apart from the computational gain, the main advantage of 
using Eq.[3]is to define a proportional relation between the clo- 
sure phase and the companion contrast. This reduces the number 
of parameters to be explored and allows the results to be pre- 
sented in a synthetic way. Additionally, we note that in Eq.[3] the 
magnification factor m takes a maximum value of about 2.6 rad 
(» 149 deg) for purely geometrical reasons. This magnification 
is only achieved in the optimal geometrical case, where all in- 
terferometric baselines add together to amplify the companion 
signal. Therefore, as a quantitative example, we can already con- 
clude that reaching a dynamic range of 10~ 3 with one single clo- 
sure phase measurement requires a closure phase accuracy bet- 
ter than 0.15 deg. In practice, this accuracy is generally reached 
after appropriate uv -plane or temporal averaging of a series of 
individual closure phase measurements. 

2.2. The close-companion limit 

In the specific case where the binary is not resolved by any in- 
dividual interferometric baselines (all a < 1), it is possible to 
simplify Eq.[3]by expanding the sines for small values of a 



a 2 n a 23 + «12 aL . 
+ o(a ) 



(4) 



We can arbitrarily choose an to be associated with the most 
resolving baseline and 0-23 with the least resolving baseline in 
the direction of the considered binary. We introduce the r ratio 
of the two geometrical terms: Q-23 = -ran (the negative sign 
is introduced for the triangle to close with positive values of r). 
Rewriting Eq. |4] we obtain 

4> 1 3 

P L 

This equation already shows two interesting aspects: (i) the clo- 
sure phase signature of an unresolved companion is proportional 
to the baseline length at the third power, and (ii) the arrays that 
are redundant when projected onto the binary direction (r = 0.5) 
seem optimal for detecting a close companion using a given long 
baseline. 



2.3. Deriving sensitivity limits 

A typical interferometric observation is not composed of a sin- 
gle closure phase measurement. We now derive sensitivity limits 
in terms of companion contrast for a set of n linearly indepen- 
dent closure phase measurements <p, of accuracy cr,. We base 
our analysis on the x 2 of the data with respect to a model of 
an unresolved source with no companion, for which all closure 
phases are zero 



X' 



Li a2 ' 



(6) 



i=i 



1 Using m apertures, one can form (m - \)(m — 2)/2 linearly inde- 
pendent closure phases. This is equivalent to holding one aperture fix ed 
and forming all possible triangles with that aperture (Monnier 2003). 
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We introduce the simplification presented in Eq. [3] and assume 
that all the measurements have similar accuracies (a valid as- 
sumption when observing an unresolved target) to obtain 



X 



z " 

p - y 



(7) 



The probability that the data set is compatible with the single- 
star model is given by the complement of the cumulative proba- 
bility distribution function (CDF) with n degrees of freedom 



P = 1 - CDF„et 2 ) 



(8) 



If P is below a predefined threshold, the data set allows the 
model with no companion to be rejected. The threshold can gen- 
erally be fixed at a 3-cr level, i.e., at a probability of 0.27%. 
The choice of the threshold actually depends on the context of 
the observations. If a large region of parameter space (or in- 
deed number of targets) is searched for a companion, then a 5- 
cr threshold may be more appropriate. A higher threshold may 
even be needed in the case of sparse data sets dominated by non- 
Gaussian systematic errors. Taking the 3-cr level as an example, 
this threshold can be converted into a sensitivity limit in terms 
of companion contrast 



CDF:'(1 -0.27%) 



(9) 



where CDF" (1 - 0.27%) is simply the \ 2 value for a 3-cr de- 
tection with n degrees of freedom. As a consequence of the ap- 
proximation presented in Eq. [3] the sensitivity limit is directly 
proportional to the accuracy of the closure phase measurements. 

Because of the sparse structure of the point spread func- 
tion associated with the diluted aperture of an interferometer, the 
depth to which a companion can be detected strongly depends on 
the relative orientation between the companion and the interfer- 
ometric baselines (information embedded in the magnification 
factors mi). For some lucky separations, the greatest dynamic 
range is achieved, while in the worst cases even obvious binaries 
with equal brightnesses can be missed. Consequently, the sensi- 
tivity limit should be defined as a two-dimensional map or for 
various completeness levels on a given search region. 



50- 



-50- 




Fig. 1. Interferometric configuration offered at VLTI. The con- 
figurations using the four relocatable Auxiliary Telescopes are 
represented by colours. The configuration using the four fixed 
Unit Telescopes is represented in black. 



Finally, Eq. Q] assumes that the angular diameters of the bi- 
nary components are unresolved by the interferometric base- 
lines. Resolving the diameter of the faint component indeed ap- 
pears unrealistic, although maybe not for the central star, es- 
pecially in the case of bright late-type giants or very nearby 
stars. In this latter situation, Eq. [TJ underestimates the closure 
phase signal, which peaks for a fully resolved primary star. This 
effect, referred to as closure phase nulling, can lead to larger 
magnification factors than the limit m < 149 deg presented in 
S ect. 12. II This specif ic observing technique is discussed in detail 
inlChelli et al.| (l2009t) . while on-sky a pplica tions can be found in 
iMonnier et all (120061). [Lacour et all d2008l) . IZhao et al.l (120081) . 
and lDuvert et al.1 d20101) ~ 



2.4. Validity limit of this study 

In the case of wide companions, care should be taken regarding 
the chromaticity limit of our study. The results presented here are 
indeed formally valid only for monochromatic light. The main 
effect of wavelength smearing inside Eq.fljis to degrade the dy- 
namic range. To avoid significant smearing, the spectral resolu- 
tion should be higher than the a quantities defined in Eq.|2] This 
translates into a spectral resolution R > 10 for a 40 mas binary 
observed with a 100 m baseline at 1.7 fim (H band). Such a low 
spectral resolution is available in most modern interferometric 
beam combiners. In the case of spatially filtered beam combin- 
ers, a similar effect may occur because of baseline smearing, in 
the case where the telescope size cannot be neglected in Eq. [TJ 
For the 1.8-m Auxiliary Telescopes, this corresponds to angular 
separations of about 175 mas for H-band observations. For 10-m 
class telescopes, this corresponds to angular separations of about 
45 mas. In practice, these limitations are not very relevant to our 
study because AO-assisted spare aperture-masking imaging on 
10-m class telescopes becomes more efficient than long -baseline 
interferometry for separations larger than about 40 mas (see e.g. 
iKraus et al.ll2008l:lLacour et al. 11201 ll) . 



3. Optical interferometric array 

We use the formalism introduced previously to compute and 
compare the capabilities of various four-telescope interferomet- 
ric configurations. 

3.1. Performances of VLTI configurations 

We compute the map of the 3-cr sensitivity limit using all the 
VLTI configurations displayed in Fig. [TJ assuming a target at 
declination -35 deg. The limits are computed considering data 
sets consisting of respectively one pointing at an hour angle 
HA = -2 h, three pointings at hour angles HA = -2 h, - 1 h 
and h, and five pointings at hour angles HA = -2 h, - 1 h, 
Oh, +II1, and +2h. A closure phase accuracy of 0.25 deg is 
assumed (see Sect. We consider a maximum binary separa- 
tion of 40 mas, which is the separation where AO-assisted spare 
aperture-masking imaging on 10-m class telescopes becomes 
more efficient than long-baseline interferometry. 

Detailed results of the widest AT configuration (A1-K0-G1- 
II) are displayed in Fig.[2]for illustration. The figure shows that 
the detection completeness for a given sensitivity level increases 
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Fig. 2. Leff: Map of the 3-cr sensitivity with the A1-K0-G1-I1 configuration from VLTI. Radial zones are the bins in separation used 
for the plots in the right panel. Right: Sensitivity as a function of angular distance, for various completeness levels. From top to 
botom: Simulations for a single snapshot pointing (top), and for three pointing (middle), and for five pointing (bottom). The contrast 
axes can be scaled for any accuracy on the closure phase (here <x = 0.25 deg). 



drastically with the number of pointings. This is clearly illus- 
trated by the decreasing number of "blind spots" (white zones) 
in the left-hand side plots. The sensitivity level is mostly flat 
for angular separations larger than 2 mas, while the detection 
performance drops considerably within this inner working an- 



gle (IWA). The median sensitivity levels in the region 2-40 mas 
are respectively about 6 x 10~ 3 , 4.5 x 10~ 3 and 4 x 10~ 3 for 
the three considered data sets. For more than five pointings, the 
median sensitivity level would continue to improve slightly, but 
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Fig. 3. Left: 3-cr sensitivity versus separation for a completeness level of 80%. Right: Completeness in the separation range 6-40 mas 
versus contrast. Simulations are for a single snapshot pointing (dashed lines) and for five pointings separated by one hour (solid 
lines). Colours are for the four configurations of VLTI displayed in Fig.Q] The contrast axes can be scaled for any given accuracy 
of the closure phase measurements (here cr = 0.25 deg). 



the shape of the sensitivity curve would no longer significantly 
change. 

The relative performances of the different configurations are 
illustrated in Fig. [3] (left: sensitivity versus separation, right: 
completeness level versus contrast). All configurations provide 
flat performances for large separations, down to their respective 
IWA where the performances dramatically drop. The IWA are 
respectively 2 mas for the A1-K0-G1-I1 and U1-U2-U3-U4 con- 
figurations, 3 mas for the D0-H0-I1-G1 configuration, and 6 mas 
for the A1-B2-C1-D0 configuration. They correspond to the spa- 
tial resolution of the smallest baseline of the array. Given the 
similarity of the results for all configurations, we discuss them 
together in two different regimes: (i) the close-companion and 
(ii) the wide-companion regimes. 

3.2. Close-companion regime 

Close companions are defined here as companions with angular 
separations that are not fully resolved by at least one baseline 
(that is B ■ \/A < 1). In this regime, the achievable contrast for 
a given completeness follows a power law of the angular sepa- 
ration p oc A~ 3 , as predicted by Eq.0 The exact factor entering 
into this law depends on the array geometry but, as expected, the 
longest arrays provide the highest spatial resolution and are thus 
able to detect both the deepest and the closest binaries. 

This well-known result was discussed by Lacha umd (l2003h 
in the context of partially resolved interferometric observations. 
We emphasize that our study additionally provides a quantita- 
tive estimation. As a typical example, we now detail the case 
of a faint companion with a contrast of 5 x 10~ 3 . We first con- 
sider the companion to be located at 2 mas from the central star. 
A closure-phase accuracy of 0.25 deg results in a detection effi- 
ciency of about 50% using three pointings with the configuration 
A0-K0-G1-I1, according to Fig. [2] (middle-right plot). However, 
if we now consider the companion to be located at 1 mas from 
the central star, the closure-phase accuracy should be 0.025 deg 
to reach the same efficiency. Since the angular separation is only 
marginally resolved by the interferometer, the lack of spatial res- 



olution has to be compensated for by an increase in the accuracy 
on the signal (super-resolution effect). 



3.3. Wide-companion regime 

We note that those companions are wide only in the interferomet- 
ric sense, corresponding to separations larger than about 4 mas 
for the typical ~ 100 m baselines available in modern interfero- 
metric facilities. 

In this regime, the detection efficiency becomes independent 
of the companion separation. Interestingly, all arrays have the 
same efficiency. In other words, as long as the companion is ex- 
pected to be resolved by the interferometric baselines, the choice 
of array configuration does not matter. We conclude that there 
is no reason to favour a given VLTI configuration when look- 
ing for faint unknown companion with separations in the range 
6 - 40 mas. More quantitatively, Fig. [3] (right) displays the de- 
tection efficiency in this annular region for the four VLTI con- 
figurations versus the companion contrast, and for two observ- 
ing scenarios (snapshot and long integration). The combination 
of five observations separated by one hour provides a detection 
efficiency higher than 95% for companion contrasts of 10~ 2 , as- 
suming a realistic closure phase accuracy of 0.25 deg. 

We note that the curve of completeness versus contrast be- 
come significantly sharper when accumulating observations. As 
shown by the solid lines in the right panel of Fig. [3] when ac- 
cumulating five pointings, the efficiency drops from 80% for a 
contrast of 5 x 10~ 3 to less than 10% for a contrast of 3 x 10~ 3 . 
The constraints provided by this dataset can thus be presented as 
a sensitivity limit and an inner working-angle, as for a classical 
imaging observation. 

Quantitatively, when accumulating several pointings, these 
detection limits computed from the derivation of Sect. [2] are 
compatible with the blind-test analyses presented by Absi l et al. 
(2011, Fig. 5) and the Monte-Carlo simulations of lLacour et al. 
00 111 Fig. 4 and Eq. 2). 
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Fig. 4. Fake VLTI interferometric configurations used in this pa- 
per: non-redundant linear D0-E0-H0-K0 (red), redundant linear 
D9-G2-H9-K9 (green, fake stations), and Y-shaped E0-J3-J2-H0 
(blue). 

3.4. Performances of standard configurations 

To add generality to the results presented in the previous section, 
we now study configurations that are not specifically linked to 
any existing interferometric array. We select the configurations 
presented in Fig. |4] which is a non-redundant linear configura- 
tion, a fully redundant linear configuration and a Y-shaped con- 
figuration. All configurations have their longest baseline of the 
same size. The results are the following: 

1 . The linear non-redundant and Y-shaped configurations have 
similar detection limits as the currently offered (irregular) 
VLTI configurations for snapshot observations. Surprisingly, 
they also have similar performances when accumulating sev- 
eral pointings, while we may have expected that the Y-shape 
would unveil faster the remaining blind-spots. 

2. For snapshot observations, the linear redundant configura- 
tion favours the dynamic range with respect to the complete- 
ness: it has a fainter detection limit for completeness lev- 
els below 50%, but becomes significantly worse for higher 
completeness levels. When accumulating several pointings, 
both the highest completeness and largest dynamic range are 
reached, although the gain is never higher than 20%. 

3. Y-shaped arrays have smaller inner working angles than lin- 
ear configurations of identical maximum baseline length, 
even when considering the accumulation of several point- 
ings. The gain is almost a factor of two in terms of angular 
resolution. 



4. Closure phase accuracy and achievable dynamic 
range 

4. 1 . Photon and piston noises 

We consider the theoretical photon noise limit for an observa- 
tion of 1 h on a star of sixth magnitude using one-metre class 
telescopes (e.g., the 1.8 m auxiliary telescopes of the VLTI). The 



choice of a sixth-magnitude star is driven by the current sensi- 
tivity limit of most interferometric instruments world-wide. A 
crude estimation of the photon noise is given by 

360 deg 

CTphot ~ =— , (10) 

V# 

where N is the total number of detected photons. If we consider 
that the 1 h observing time should include the overheads and the 
observation of a calibration star, the effective integration time on 
target will be of the order of 20 min. Assuming a realistic total 
transmission of 2%, including both the instrumental and atmo- 
spheric contributions, the number of detected photons is N ^ 10 s 
and the resulting photon noise is cr ; ,/, 0( ~ 0.1 deg. For a l mag star, 
the resulting photon noise would be <x ; ,/, 0( ss 0.01 deg. 

In theory, closure phase is a robu st observable against the 
telescope phase errors dMonn ier 2003). However, in the context 
of non-zero exposure times and the presence of atmospheric tur- 
bulence, closure phase measurements are also affected by piston 
noise. A proper estimation of its amplitude is beyond the scope 
of this paper but it is still possible to provide a rough upper limit. 
Since piston noise is independent of the number of incident pho- 
tons, it is expected to dominate the final statistical uncertainty 
for very bright stars. As an example, in the case of the PIONIER 
instrument at VLTI, the statistical uncertainty for bright stars is 
typically of the order of 0.25 deg to 2.5 deg for an integration 
time of 1 min. This uncertainty depends on the atmospheric con- 
ditions as expected for piston noise. In decent atmospheric con- 
ditions, integrating over 20 min allows the piston noise contribu- 
tion to be reduced below 0.2 deg, as it decreases with the square 
root of the integration time. 

4.2. Calibration accuracy 

It is interesting to compare these fundamental limits to published 
accuracies that include the calibration of the instrumental closure 
phase (also called the transfer function): 

VLTI/ AMBER: lAbsil et alJ (1201 Oh reported calibration errors of 
between 0.20 deg and 0.37 deg depending on the night, us- 
ing this three-telescope combiner in its medium spectral res- 
olution mode (R = 1 500). With the low spectral resolu- 
tion mode (R = 35), typical calibration errors range from 
one to a few degrees (see for instance iKraus et aL 2009; 
iLe Bouquin et alJl2009h . 

VLTI/PIONIE R: Typical calibration errors range from 0.25 deg 
to 1 deg dLe Bouquin et alj|201ll lAbsil et all 1201 lb for this 
four-telescope combiner. Sequences with closure phases sta- 
ble down to 0. 1 deg have been recorded. Systematic discrep- 
ancies have been noted when calibration stars were separated 
by more than 10 deg on the sky. 

CHARA/MIRC: The typical accuracy obtained with this four- 
telescope combiner is between 0.1 and 0.2 deg, which makes 
this instrument the most accurate of the currently available 
suite. Calibr ation uncertainties dominate the final accuracy 
at this level (IZhao et alJ[2008L[20ToL[20in) . 

Altogether, a typical noise floor of ~ 0.25 deg seems to 
appear for the calibration of the closure phases in long base- 
line interferometry. Two results indicate that the major cause 
is probably longitudinal dispersion: (i) that the accuracy de- 
pends on the spectral resolution and (ii) the dependence on po- 
sition of the calib r ation star on the sky. This is also the find- 
ing of IZhao et alJ (1201 ll) . who proposed an elaborate calibra- 
tion scheme for MIRC. Although this is clearly a very promis- 
ing way of characterizing already known substellar companions, 
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Fig. 5. Left: 3-cr sensitivity versus separation for a completeness level of 80%. Right: Completeness in the separation range 6-40 mas 
versus contrast. Simulations are for a single snapshot pointing (dashed lines) and for five pointing separated by one hour (solid lines). 
Colours are for the four configurations displayed in Fig. |4] The black curves are for the U1-U2-U3-U4 configuration displayed in 
Fig.Q] The contrast axes can be scaled for any accuracy on the closure phase (here cr = 0.25 deg). 



this method is probably not suited to surveying a large num- 
ber of stars with a standard calibration proc edure. Interestingly , 
0.25 deg is also the noise floor reported by lLacour et al. (2011) 
for the calibration of the closure phase of the spare aperture 
masking mode of NACO at VLT. This calibration noise floor 
of 0.25 deg theoretically does not prevent us from reaching very 
high dynamic ranges, by accumulating a large number of obser- 
vations and/or baselines, as for instance in sparse aperture mask- 
ing or spectrally dispersed observations. Care should however be 
taken to ensure that individual closure phase measurements are 
statistically independent. In particular, one should avoid repeat- 
ing the same systematic errors in individual data sets, e.g., by 
choosing different calibrator stars, instrumental setups, etc., in 
order not to reach a true noise floor in the observations. 

Concerning future instruments, the announced accuracy on 
the closure phases is 1 deg for the K-band four-telescope com- 
biner GRAVITY (d ocument VLT-SPE-ESO-15880-4853 and 
iGillessen etai1l2010l) and from 1 deg to 5 deg for the L-band 
four-telescope com biner MATISSE (F lorentin Millour, private 
communication and lLopez et al.ll2008l) . Although these perfor- 
mances may be conservative, we conclude that the next genera- 
tion of VLTI instruments is unlikely to break the 0.25 deg limit. 

4.3. Discussion 

It appears realistic to reach a closure phase accuracy of 0.25 deg 
within less than one hour on one-metre class telescopes for 
stars of magnitude six and brighter. According to the results of 
Sect. [3] such performances allow a dynamic range of 5 x 10~ 3 
(Amag = 5.75) to be reached with 80% completeness when five 
pointings are obtained with a four-telescope interferometer. The 
same performance would probably be reached within a snapshot 
using an interferometric instrument combining six telescopes or 
more at a time. 

This result can be compared with the survey for stellar and 
sub-stellar companions on the ten-metre Keck and five-metre 
Palomar telescopes using apertu re masking techniques in the 
K band dKraus et"aT]l2008L 1201 lh . The achieved dynamic range 



is AK m 5.5 for separation as small as 25 mas (slightly worse 
for Palomar). A similar dynamic range and inner working an- 
gle have been achieved within an ongoing survey of massive 
stars using aperture masking at VLT/NACO in the H band, e.g. 
AH « 5 down to 25 mas with this eight-metre diameter telescope 
(Hugues Sana, private communication). All close companions 
presented in these near-infrared surveys would have been de- 
tected by interferometry with an efficiency higher than 90%, pro- 
vided that they reside within the interferometric field-of-view. In 
addition, this efficiency would have been achieved down to about 
2 mas. At longer wavelengths, the ape rture masking techn ique 
has a dynamic range of about AL a: 7.5 (Hi nklev et al.ll201 lb . al- 
though the inner working angle in that case is only 70 mas. There 
is currently no L-band interferometric beam-combiner with clo- 
sure phase capabilities to which these performances could be 
compared. 

Several observing programs related to faint companion de- 
tection would benefit significantly from the capabilities of clo- 
sure phase measurements on long-baseline interferometric in- 
struments. An example is the search for low-mass (sub)-stellar 
companions around main-sequence stars residing in nearby 
young associations. Considering associations with ages between 
10 Myr and 200 Myr, and a limiting magnitude K = 6 for the in- 
strument, one could survey stars up to about 15-20pc for stellar 
type MOV, 40 pc for type GOV, and 120pc for type AOV With 
an estimated median dynamic range of AK 6, we computed 
the masses of the faintest companions that could be detected 
within a survey of n earby moving groups, usin g the (sub-)stellar 
cooling models of Baraf fe et al.l (Il998l 120031) . The results are 



given in TableQ] showing that the 13M Jup (= O.O12M ) limit be- 
tween the brown dwarf and planetary regimes can be reached for 
young late-type dwarfs. In the case of A-type stars, one would 
be sensitive to companions in the range M3V-M7V depending 
on the age. For even younger stars, located in nearby star form- 
ing regions, closure phase measurements have the potential to 
reveal the forma t ion of planetary-mass objects, as suggested by 
lKraus& Ireland! (l2012h . 
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Table 1. Sensitivity limits in terms of companion masses around 
young main-sequence stars. 



Age 


AOV 


GOV 


MOV 


lOMyr 


0.09M o 


0.017M o 


0.012M o 


50Myr 


0.22M o 


0.043M o 


0.013M o 


200 Myr 


0.35M o 


0.08M o 


0.030M o 



Another program is the determination of the binary fraction 
for massive stars. The interest is that despite the preponderance 
of multiple stars, the mechanism that produces multiple stars 
rather than single stars is still uncertain. The measurement of the 
mass distribution, and how it evolves with the mass of the pri- 
mary, is an appropriate tool for disentangling between capture 
and fragmentation models. Stellar co mpanions to B-type stars 
have been investigated using AO (e.g. iRoberts et alJl2007l) . al- 
though the stars observed typically have a large range of dis- 
tances, limiting the statistical significance of the results. Radial 
velocity measurement of massive stars is challenging owing to 
the lack of suitable spectral lines and their intrinsic broadening. 
With the limiting magnitude K — 6 presented in this paper, it 
is possible to observe interferometrically all the B stars within 
a distance of 75 pc (~100 objects for the southern hemisphere), 
providing the first comprehensive study of massive binaries in 
the 0.25 - 5 AU separation range. 

Last but not least, one of the main selling argument for high- 
precision closure phases in optical interferometry has been the 
direct detection of hot extrasolar giant planets (EGP). Several 
hot EGP host stars are indeed bright enough to be observed with 
state-of-the-art interferometric instruments. For mature plane- 
tary systems (> 10 Myr), the expected contrast between the 
planet and the star is however generally too low (< 10~ 3 ) 
to be currently accessi ble with closure phase measurements 
(IZhaoetal .1120081 1201 lh . To routinely reach the hot EGP regime 
(AK =i 8- 10), a gain of two to four magnitudes is required in the 
dynamic range, which would translate into a noise floor of be- 
tween 0.04 deg and 0.006 deg on the closure phase. Achieving 
such an accuracy would probably require a significant break- 
through in the instrumental domain. 

5. Conclusions 

In summary, optical interferometric surveys designed to detect 
faint companions have the following properties: 

1. The observable (closure phase ) is robust ag ainst unstable 
atmospheric seeing conditions (Monnier 2003). Integrating 
over 20min is sufficient to consistently reduce the photon 
and atmospheric noises below 0.25 deg, which appears as a 
hard limit for the calibration of current instruments. 

2. A single snapshot with four telescopes provides a 80% de- 
tection efficiency at Amag = 4.5 as soon as the binary sep- 
aration is fully resolved. The only requirement of the inter- 
ferometric array is to use baselines as long as possible to 
improve the inner working angle, which is typically of the 
order of a few milli-arcseconds. 

3. Accumulating more observations (several pointing and/or 
recombining more telescopes) allows a dynamic range 
Amag = 6 to be reached, which appears to be a realistic limit 
in respect to published performances. Going deeper would 
require us to break the current limit of 0.25 deg on the clo- 
sure phase accuracy, or to massively increase the number of 
observations. 



4. The achievable dynamic range scales linearly with the clo- 
sure phase accuracy. 

In conclusion, interferometric closure phase surveys would 
be well-suited as filler programs for service-mode interfero- 
metric facilities, such as the VLTI. They can be considered as 
a useful complement to the AO-assisted imaging surveys cur- 
rently carried out on ten-metre class telescopes. In particular, 
the search space of long-baseline interferometry bridges the gap 
between the wide companions found in direct imaging and the 
close companions detected by RV measurements. Moreover, in- 
terferometry could nicely complement RV studies in the par- 
ticular cases where RV measurements are quite inappropriate. 
Young stars, for instance, are especially promising targets since 
their (sub)stellar companions are supposed to be relatively bright 
compared to their host stars. 
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